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r> ! ABSTRACT 

D ' Semi-analytic treatment of the power spectrum with the approximation of constant 

|JLh I hnear bias provides a way to compare cosmological models to a large amount of data, as 

P^ , Peacock & Dodds (1994, 1996; hereafter PD94 and PD96, respectively) have shown. 

^__l ' By applying the appropriate corrections to the observational power spectrum it is 

possible to recover the underlying linear power spectrum for any given cosmological 

p<| ■ model. Using extensive N-body simulations we carefully test and calibrate all impor- 

K*" ' tant corrections. To demonstrate that the method is applicable to a wide range of 

On , cosmological models, we tested the standard D, — 1 cold dark matter (CDM) model, 

QN ■ O < 1 models that include a cosmological constant (ACDM), and Q = 1 models with 

a mixture of cold and hot dark matter (CHDM), both with one massive neutrino and 
two equal mass neutrinos. The ACDM and CHDM models are normalized to COBE 
and to cluster abundances. Our tests indicate that the improved linear-nonlinear map- 
Qv I ping recently suggested by PD96 for treating CDM-type power spectra works well for 

■~>^ I a wide range of scale-dependent power spectra. However, we find that the recovery 

i-C ' of the linear power spectrum from observations following PD94, which is often used 

'T^, to test cosmological models, can be misleading because the corrections are model- 

Q ■ dependent. A model should not be judged based on the linear spectrum recovered by 

}—{ I that procedure, but must be compared with the recovered spectrum for that particular 

C/3 , model. When we apply the proper corrections for a given model to the observational 

Cu ' power spectrum, we find that no model in our test group recovers the linear power 

^ I spectrum well for the bias values suggested by PD94 between Abell, radio, optical, 

and IRAS catalogs: bA '■ b-^ : bo ■ h — 4.5 : 1.9 : 1.3 : 1.0, with bi = 1.0. The recovered 
linear ACDM and CHDM power spectra were systematically below their respective 



o 

O 



H , linear power spectra using bi = 1.0. When we allow 6i to vary (keeping the same bias 

ratios) we find that: (i) CHDM models give very good fits to observations if optically- 
selected galaxies are slightly biased {bo ~ 1.1). (ii) Most ACDM models give worse 
but acceptable fits if blue galaxies are considerably antibiased (0.6 ^ bo ^ 0.9) and 
fail if optical galaxies are biased, (iii) There is a universal shape of the recovered linear 
power spectrum of all ACDM models over their entire range of explored wavenumbers, 
0.01 ^ k ^ 0.6 h Mpc~^. For a given bias, recovered linear power spectra of CDM 
and CHDM models are nearly the same as that of ACDM in the region 0.01 ^ k f^ 0.2 
h Mpc~^ but diverge from this spectrum at higher k. We tabulate the recovered linear 
spectra, and also the initial linear spectra, for all the models considered. 
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1 INTRODUCTION 

The power spectrum of density fluctuations P{k) is a pow- 
* E-mail: erful tool for investigating the large-scale structure of the 

aklypin@nmsu.edu (AK); gross@physics.ucsc.edu (MAKG); Universe. It is a useful statistical test of the distribution 
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larger than superclusters, and it discriminates between cos- 
mological models (see, for example, Liddle et al. 1996ab; 
Ma 1996; Klypin, Primack, & Holtzman 1996, hereafter 
KPH96). However, there are several obstacles that must be 
overcome to relate the observed nonlinear power spectrum to 
theoretical linear power spectra. Because observations typ- 
ically give estimates for luminous objects in redshift space 
while theoretical models describe dark matter in real space, 
the differences between luminous and dark matter must be 
taken into account before a direct comparison between data 
and theory can be made. 

Since the radial positions of galaxies are determined 
from their observed recession velocities, any peculiar mo- 
tions (from either linear or nonlinear growth of fluctuations) 
will distort their placement. Kaiser (1987) estimated the 
effect for long waves in the linear regime and included a 
scale-independent bias. This correction applies to the region 
k ^ 0.1 h Mpc~^ but not to smaller scales which are in the 
nonlinear regime. To treat this region, Peacock (1992) used a 
simplified model to estimate the effect of random motions in 
collapsed objects on the observed power spectrum. Peacock 
& Dodds (1994, hereafter PD94) combined the two separate 
effects, along with the overall bias level, into a single conve- 
nient expression that is claimed to remove these effects from 
the nonlinear power spectrum. 

Additional nonlinear effects must be accounted for as 
well. Theory provides a linear power spectrum for a cos- 
mological model, but the spectrum of the real universe has 
undergone significant modification because of nonlinear evo- 
lution. Numerical simulations are useful for estimating the 
nonlinear power spectrum but are time consuming and com- 
putationally expensive. A more practical approach would be 
to devise a general method for analytically constructing the 
nonlinear power spectrum from the linear power spectrum. 
This would allow any model to be evaluated quickly and 
open the possibility of directly removing nonlinear effects 
from a real power spectrum. Throughout this paper we will 
refer to the method of recovering the linear spectrum from 
the raw observational spectrum by applying a set of cor- 
rections as the 'linearization' procedure. We will refer to the 
inverse of linearization, the method of predicting the nonlin- 
ear spectrum of galaxies starting from the linear dark matter 
spectrum, as the 'nonlinearization' procedure. 

Several advances in developing linearization methods 
have been made in the last few years. By considering the 
collapse of spherical structures and the nonlinear evolution 
in the limit of stable clustering, Hamilton et al. (1991, here- 
after HKLM) devised an analytical method for relating the 
linear correlation function on a linear scale to the nonlinear 
correlations on a smaller, nonlinear scale. PD94 extended 
this method to power spectra, allowed for open models and 
models with a cosmological constant, and provided a fitting 
formula optimized by a series of numerical simulations. This 
semi-analytical method worked reasonably well for scale-free 
power-law spectra, but failed for models with spectral in- 
dices n ^ —2 including scale-dependent spectra like those 
of CDM models that tend to n ~ —3 on the smallest scales. 
Like HKLM, PD94 expected that their series of corrections 
would be nearly model-independent and could be applied 
to observational power spectra to arrive at the underlying 
linear power spectrum of the Universe. Not long after, how- 
ever, Jain, Mo, & White (1995) found that there is indeed 



a significant dependence on the model through the spectral 
index n of the linear power spectrum and suggested that an 
'effective index' be used to treat curving spectra like those of 
CDM or CHDM models. This prompted Peacock & Dodds 
(1996, hereafter PD96) to revise their earlier fitting formula 
to include dependence on n, and propose that n should vary 
with scale to treat curving spectra. 

In this paper we adopt the analytical form of the red- 
shift corrections of PD94 and the improved linear-nonlinear 
semi-analytical approximation of PD96, although we find 
it necessary to modify both of these slightly. We have ex- 
tensively tested each step of the procedure using numeri- 
cal simulations of unbiased dark matter in CDM, ACDM, 
and CHDM models. We find that the redshift corrections of 
PD94 are close to those that we determine from our sim- 
ulations, although they are improved when we use the ap- 
propriate ID velocity dispersions for each model instead of 
universally applying 300 km s~^. We also demonstrate that 
our modification of the linear-nonlinear mapping of PD96 
works very well for the cosmologies that we examine; CDM, 
ACDM, and CHDM models are all well fit. 

Having found that the method can handle a variety 
of models, we applied it to a compilation of observational 
power spectra (APM, APM/Stromlo, CfA-2, QDOT, and 
radio galaxies, and Abell clusters). Our intended goal is to 
demonstrate that the linear power spectrum can be recov- 
ered from concrete catalogs of objects in redshift space or 
real space for a number of interesting cosmological models 
and then use this technique to determine which model best 
fits the observed power spectrum. Because the corrections 
are model-dependent (mostly through the shape of the linear 
power spectrum), no simple scheme can reproduce the un- 
derlying linear spectrum without detailed information about 
the model. This dictates that we must select well-defined 
models and try to reconstruct their linear power spectra 
using the appropriate corrections for each model. For this 
reason, the widely used power spectrum of PD94 is of lim- 
ited utility, since no model was specified. Several papers have 
used this power spectrum as a means for testing cosmological 
models (e.g., Liddle et al. 1996ab; Ma 1996; Coble, Dodelson, 
& Frieman 1996). We assert that this is incorrect and can 
lead to incorrect conclusions. The only measure of success 
of a particular model in this context is its ability to recon- 
struct the linear power spectrum from the model-dependent 
linearization of the observational power spectrum. 

We have followed PD94 in assuming that the bias is lin- 
ear and constant, and in our choice of observational data, 
but we have propagated the errors somewhat differently. We 
also considered a form of scale-dependent bias for one model. 
We have used the bias ratios of PD94, but we allowed the 
reference level bi to float and by adjusting the overall bias 
level we optimized the fit of the linearization to the linear 
power spectrum for each model. We find that the lineariza- 
tion of the observational power spectra according to the 
ACDM family of models can all be described by the same 
shape and amplitude for a given bias level. In contrast, lin- 
earizations of CHDM models do not have a common form. 
Moreover, it appears that there is no universal linearized 
power spectrum since CDM, ACDM, and CHDM lineariza- 
tions are all easily distinguishable. In addition, we arrive at 
the conclusion that optically selected galaxies must be an- 
tibiased {bo — 0.6 — 0.9) if they are to be compatible with 
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the ACDM family of models. The estimates of the optical 
biases of different models from potent analyses by Hudson 
et al. (1995) seem to indicate that a small amount of an- 
tibias is acceptable for some models, though Peacock (1996) 
would disagree. Optical galaxies in CHDM models must be 
slightly biased (60 — 1-1), which would probably tend to 
favor CHDM models. 

This paper is organized as follows. The details of our 
simulations are described in §B. §0 lists the various power 
spectra in the observational data set. In §M we review the 
semi-analytic method in depth. We discuss the tests of the 
method and its application to our simulations in §0 The lin- 
earized observational power spectra for all the models that 
we consider are presented in §M. Our results are summarized 
in ^. Appendix |A| gives theoretical power spectra for all the 
models considered. 



2 NUMERICAL SIMULATIONS 

In order to fully test our semi-analytical approximations, we 
performed an array of numerical simulations of several dif- 
ferent cosmological models. All simulations were performed 
using standard Particle-Mesh (PM) codes (Hockney & East- 
wood 1981; Kates, Kotok, & Klypin 1991; Smith 1995; 
Klypin, Nohhenius, & Primack 1997; Gross et al. 1997). Ta- 
ble fl lists the models and gives their defining parameters. 
We used the 'standard' cold dark matter model (SCDM) 
with bias 6 = a^^ = 1.5, four different flat ACDM mod- 
els, and two n = 1 CHDM models. The CHDM0.7 model 
has one neutrino species and the CHDMo.g model has two 
equally massive neutrino species. Except for SCDM, all our 
simulations started from COi?-E-normalized Gaussian initial 
conditions. 

Table H shows the box sizes, resolutions and number of 
realizations for each model. Multiple box sizes and realiza- 
tions were simulated for many models in order to reduce 
the impact of two effects which can change estimates of the 
nonlinear power spectrum as calculated from simulations. 
(1) Cosmic variance can cause significant fluctuation of the 
power spectrum on large scales due to poor counting statis- 
tics. Averaging the power spectra of several smaller realiza- 
tions together reduces this effect significantly. (2) Small box 
sizes (L ^ 100 h~^ Mpc) have the effect of reducing the 
power in the nonlinear region of the power spectrum by as 
much as 50 per cent, leading to an underestimation of the 
scale at which nonlinear effects begin to appear. This is due 
to the fact that long waves (which are coupled to shorter 
waves) cannot be present in the box. Boxes with L ^ 200 
h~^ Mpc are better for constructing accurate power spectra 
of the quasi-nonlinear region because the first several bins 
at small k then lie in the truly linear regime. By using boxes 
with L — 300 h~^ Mpc for most of the models discussed in 
this paper, we address both problems. 

For CDM-type spectra like those that we examine in 
this paper, (with or without a cosmological constant, but 
not including CHDM), the primordial spectrum is modi- 
fied by physical processes that can be expressed in a scale- 
dependent transfer function, T{k). Since the shape of the 
CDM spectrum does not change much after the epoch of 
equality, it is a function of k only. The power spectrum of 



the tilted CDM-type models considered here can thus be 
approximated as 



P(fe) = Afe"T^(fc) 



2.,.^5n"W) 



g2(no) 



(1) 



where g{^, SIa) is the growth rate of fluctuations used to 
specify the power spectrum at different epochs and T{k) is 
the transfer function. We use the approximation of Carroll, 
Press, & Turner (1992) (cf. also Lahav et al. 1991), 



ff(n,nA) = -n 



n*/' - Oa 



_V1 -I- ' 
2^^^ 70' 



(2) 



where f^A = A/{3Ho). The SCDM transfer function that we 
used was the commonly used approximation of Bardeen et 
al. (1986, BBKS), which has the form 

ln(l + 2.34g) 

(3) 
X 1 1 + 3.89q + {IG.lqf + {5A6qf + (6.71g)*] 



T{q) 



-1/4 



where q = k/(hV) and F is the 'shape parameter' deflned by 
Sugiyama (1995) (but cf. Hu & Sugiyama 1996 for improved 
spectra). 



f7o/iexp(— fib — f^b/^o)- 



(4) 



Here, f2b is the contribution of baryonic matter to the cos- 
mological density. However, we only used the BBKS trans- 
fer function with F = 0.5 for SCDM. For our other sim- 
ulations, the transfer functions of the ACDM and CHDM 
models we used were calculated from a full Boltzmann treat- 
ment updating that in Holtzman (1989); fitting formulas are 
given in Appendix M. For CHDM models, the different time- 
dependent growth rates of the cold and hot species means 
that the shape of the power spectrum changes over time. In- 
corporating the growth rates into the transfer function gives 
a general form of the CHDM power spectrum. 



PcHDM(fc)=ylfc"r'(fc,f7(t)). 



(5) 



The CHDM transfer functions were essentially those of 
Klypin et al. (1993) and Primack et al. (1995); cf. Ap- 
pendix \M. 

In all, we have three classes of dark matter models 
with a total of seven distinct models. Our sampling cov- 
ers a wide range of power spectra with shallow (SCDM), 
medium (ACDM), and steep (CHDM) spectral indices on 
small scales, which allows us to test our method thoroughly. 



3 OBSERVATIONAL DATA 

The observational power spectra to which we will apply our 
linearization methods consist of six independent data sets 
and two cross correlations between data sets, as assembled 
by PD94. Included in the sample are catalogs of optical clus- 
ters, radio galaxies, optical galaxies, and galaxies from the 
IRAS catalog. Each is briefly described below. 

One set of a real-space power spectrum is included, the 
APM power spectrum (Baugh & Efstathiou 1993). It is ar- 
rived at by deprojecting the angular clustering of individual 
galaxies from the APM survey. Because it only deals with an- 
gular positions, it directly yields the real-space power spec- 
trum without the need to correct for redshift-space effects. 
Of the full data set, we choose to examine only points in 
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Table 1. Summary of model parameters 



Model 


f^o" 


^a'' 


^c+b" 


n,^ 


h" 


nf 


crgS 


Q'^ 


pi 


SCDM 


1.0 


0.0 


1.0 


0.00 


0.5 


1.00 


0.67 


8.49 


0.50 


ACDMo.3 


0.3 


0.7 


0.3 


0.026 


0.7 


1.00 


1.10 


22.0 


0.20 


ACDMo.4 


0.4 


0.6 


0.4 


0.035 


0.6 


0.975 


1.00 


21.8 


0.23 


TACDMo.4 


0.4 


0.6 


0.4 


0.035 


0.6 


0.90 


0.873 


24.7 


0.23 


ACDMo.5 


0.5 


0.5 


0.5 


0.035 


0.6 


1.00 


0.818 


13.93 


0.29 


CHDMo.7 


1.0 


0.0 


0.7 


0.10 


0.5 


1.00 


0.676 


17.0 


0.39 


CHDMo.8 


1.0 


0.0 


0.8 


0.10 


0.5 


1.00 


0.719 


18.44 


0.34 



° Matter density parameter. 

Cosmological constant density parameter Ha = ^/(Sii^o)- 
'^ Cold dark matter density parameter + baryon matter density parameter. 
** Baryonic matter density parameter. 

"^ The Hubble parameter specified as h = Hq/{WO km s~^ Mpc~^). 
' Index of the power spectrum on large scales. 

^ Linear mass fluctuation amplitude (AM/M)rnis smoothed with a top-hat filter of 
radius & h~^ Mpc. 

'' Quadrupole of CMBR anisotropy QrmS-ps- 

' Sugiyama (1995) shape parameter F = f2o/iexp{— n^ — Q\,/Qo). The actual spectral 
shape of the CHDM models is steeper for large k than the ACDM models, contrary 
to the impression that the high values of F for these models might give. Note that 
the value of the modified shape parameter F = 0.5(0.95cT8/3.5ct25)^'^'^) introduced in 
Borgani et al. (1997) is nearly identical to the Sugiyama shape parameter for all but 
the two CHDM models, for which it would be 0.16 and 0.19, respectively. 



Table 2. Summary of simulations used in this paper 

Model mesh" particles'' box size'^ resolution'* realizations'^ 

size (/i^i Mpc) (/i~i kpc) 



SCDM 


768^ 


3843 


300 


390 


1 


ACDMo.3 


5123 


2563 


200 


390 


3 


ACDMo.4 


512^ 


2563 


200 


390 


3 


TACDMo.4 


768^ 


3843 


300 


390 


1 


TACDMo.4 


5123 


2563 


100 


195 


3 


TACDMo.4 


5123 


2563 


50 


98 


3 


ACDMo.5 


7683 


3843 


300 


390 


1 


CHDMo.7 


7683 


3 X 2563 


255 


332 


1 


CHDMo.8 


7683 


3 X 3843 


300 


390 


1 



'^ Total number of PM divisions. 

^ Number of particles in the box. 

'^ Total size of the box. 

'^ Box size divided by the number of PM divisions. 

'^ Number of simulations with different realizations that were averaged. 



the region 0.015 < k < 1 h Mpc~^. Power spectra from 
three catalogs of galaxies in redshift-space were used: the 
Stromlo/APM survey (Loveday et al. 1992) ; the CfA survey 
(Vogeley et al. 1992); and IRAS galaxies from the QDOT 
sample (Feldman et al. 1994). A straight mean of the two 
power spectra in the CfA paper was taken. Two catalogs of 
clusters of galaxies in redshift-space were used, the power 
spectrum of Abell clusters of richness class 1 or higher (Pea- 
cock & West 1992) and the power spectrum of radio galaxies 
from (Peacock & Nocholson 1991) The two cross-correlations 
are between IRAS galaxies and Abell clusters and between 
IRAS galaxies and radio galaxies (Mo et al. 1993). 

Different biases should be assigned to the different data 
sets based on how they were selected. PD94 defined four 
adjustable bias parameters for the four categories in this 
sample: &a for Abell clusters, bn for radio galaxies, bo for 



optically-selected galaxies, and fei for IRAS galaxies. Using 
a likelihood method, PD94 chose the 'best' ratio of these 
bias parameters to be 



feA : bn -.bo -.bi =4.5 : 1.9: 1.3: 1.0, 



(6) 



normalized to 6i = 1. In this paper we adopt the above ratio 
but use 6i as a reference level which can be adjusted. The 
necessity of this step will be discussed in §H. 

This sample well represents the observable power spec- 
trum in the quasi-nonlinear regime up to such small k that 
there should be no apparent nonlinear effects. The method 
that is described in the following section performs very well 
in this region of the power spectrum. Using this method we 
will attempt to remove all of the observational effects from 
these power spectra. 
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4 METHOD 

The linearization method that we use is whoUy based on 
that of PD94 and the later modification by PD96 which 
improved the handling of CDM-type power spectra. In or- 
der to recover a pure linear spectrum from real observations 
three effects must be corrected for: redshift distortions, bias 
and nonlinear growth. In practice, corrections for the red- 
shift distortions and bias are inseparable and are applied 
together through a single equation (§4.1). The resulting non- 
linear spectrum can be mapped to the linear spectrum on 
the linear scale. 

The contribution to the fractional density variance per 
bin of Infe, denoted A^(fe) (Peebles 1980), is related to the 
power spectrum P{k) = (|5k| ) as 



A'W = ^Pik)- 



(7) 



4.1 Correcting for redshift distortions and bias 

Kaiser (1987) gave an analytical form for the enhancement 
of the power spectrum due to the collapse of waves in the 
linear regime. For a single wave, the relation between real 
space A^ and redshift space A^ is 



A^(fc) = A?(fc) 



l + ^J.' 



f 



where / is the linear growth rate of velocity 
dlnS 



(8) 



(9) 



dm a 

/i is the direction cosine between the wave vector and the line 
of sight, and a scale-independent bias b is included. Averaged 
uniformly over all /i, this relation becomes 



A^(fc) 



A?(fc) 



, 2/ if 
1 H — — H — ^^ 
3 6 5 62 



(10) 



This formula describes the modification of the power 
spectrum on linear scales. On nonlinear scales, large random 
peculiar velocities within collapsed objects cause an appar- 
ent elongation in a direction along the line of sight ('fingers- 
of-God'), which decreases power on cluster-sized scales and 
smaller. Assuming that the velocity distribution is Gaussian, 
this modifies the power spectrum as 



A^(fc) = A,(fc)exp 



"Ho" 



(11) 



where a is the ID rms velocity dispersion (Peacock 1992). 
When averaged over all angles, this gives 



A.'(fe) 



A?(fc) 



/¥ erf (fccr/Ho) 



PD94 combined Eq. (hof) and ( |l2| ) along with the galaxy 
bias into a single formula that is applicable on all scales and 
also allows for cross-correlation between power spectra from 
two different catalogues. This formula is 

A^k) ^ Al{k)[bib2Giy,ai,a2)]-' 



where 



100 



VK+^I)72, 



(13) 



(14) 
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Figure 1. The redshift correction is plotted as a function of k 
and of the ID velocity dispersion, ct, for the TACDMo.4 model 
with 6 = 1. On large scales the overall correction is that of Kaiser 
(1987) while on small scales corrections are due to the 'finger- 
of-God' effect, which is sensitive to a. Notice that in the region 
fc <; 0.1 h Mpc~^, a slightly higher or lower value could alter the 
correction noticeably. 



^ = fm/b, 

G(y,Ql,Q2) = 

^^ [3ma2 + 2(ai + 02)2/^ + 4j/4] 



(15) 



(16) 



exp(-i; ) 



[aia2{3 + 2y^) + 2{ai + a2)y^] 



and subscripts 1 and 2 refer to quantities for two different 
catalogues, in the case of a cross-correlation. E g. (|l4| ) is for 
k in units of h Mpc~^ and b is the bias. Eq. (|l3) is to be 
applied point by point to the raw observational power spec- 
trum in redshift space. In the case of the APM spectrum, 
G must be taken as unity since there are no redshift-space 
distortions. 

Because of Eq. (h4) and (h_q), the redshift correction in 
Eq. ( |l3| ) is sensitive to the velocity dispersion. PD94 make 
the approximation that a — 300 km s~^ for each of the 
catalogs in the data set. We note that a slightly different 
choice of a could significantly affect the redshift correction 
at fc ^ 0.1 h Mpc~^. In Figure n^ we show how Eq. (nsl) ap- 
pears when several different values of the velocity dispersion 
are chosen for the TACDM0.4 model. For example, the cor- 
rection a,t k — 1 h Mpc~^is about 40 per cent greater for 
a = 400 km s'^ than for a = 300 km s"^ 

We are able to test the accuracy of the redshift correc- 
tions using our simulations. We placed dark matter parti- 
cles from the TACDM0.4 model at the z — moment into 
redshift space according to Hubble's law and the peculiar 
motions of the particles. Then we calculated the power spec- 
trum from two different vantage points within the box and 
averaged together all of these power spectra from each of 
the realizations of a particular box size. By taking the ra- 
tio of the redshift space power spectrum to the real space 
power spectrum, we produced a redshift correction for the 
model which can be directly compared with the prediction 
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Figure 2. The redshift correction vs. A; as verified with simula- 
tions. Two models were tested, (a) TACDM0.4 and (b) ACDM0.3, 
both with b = 1. The results for different box sizes (broken lines) 
are plotted along with the analytical predictions of Eq. ( 13) (solid 
and dotted lines) for two different values of a. It can plainly be 
seen that the corrections based on a taken from the simulations 
give much better fits than the straight cr = 300 km s~^ of PD94. 



of Eq. ([13|), which is presented in Figure eI Two diflerent 
theoretical predictions are shown, the PD94 prediction with 
a = 300 km s~^ and a prediction curve using the average 
velocity dispersion from simulations of the largest boxes for 
each model. Clearly, the simulated corrections conform bet- 
ter to the prediction using the actual velocity dispersion 
than the PD94's assumed value of cr = 300 km s~^. We con- 
clude that this is an important effect which must be included 
in this type of analysis, and we will therefore apply redshift 
corrections based on the ID velocity dispersions from simu- 
lations throughout this paper. 

In order to arrive at ID velocity dispersions that are 
appropriate for the models that we test in this paper, we 
examined all models that have 300 h~^ Mpc boxes. We cal- 
culated the rms velocity dispersion of all particles in the 
boxes and computed the redshift correction to the power 
spectrum as describe above. We found the following ID ve- 
locities to be representative of the three families of models 
considered: 610 km s'^ for SCDM, 430 km s'^ for ACDM, 
and 540 km s~^ for CHDM. These numbers are used for a 
throughout the rest of the paper. 



4.2 Correcting for nonlinear growth 

By extending the method of HLKM from the correlation 
function to the power spectrum, PD94 arrived at a fit- 
ting formula that would relate the linear to the nonlinear 
spectrum. This procedure has two parts: adjusting the am- 
plitude of the power spectrum, and transforming between 
the linear and nonlinear scales. Originally it was thought 
that the linear-nonlinear mapping would be nearly model- 
independent, with the density parameters being the only 
information about the model to enter into the equations (cf. 
PD94). However, later investigation by Jain, Mo, & White 
(1995) showed that there was a significant dependence on 
the slope of the linear power spectrum as well. We chose to 
adopt the new method of PD96 with small modifications for 
performing the linear-nonlinear mapping. 

We assume that the linear l^i^{k) is related to the non- 
linear ANL(fc) through 



AL(fcNL)=/NL(AE(fcL)) 

where /nl is a fitting function of the form 



/NL(a;) = X 



1 + BPx + [Ax] 



aP 



l + {[AxYg^{Q.)/[Vx^/^]Y 



(17) 



(18) 



This formula has an explicit dependence on Q, through the 
growth suppression factor g{Q,,h.), given in Eq. (S), and 
five free parameters. The best-fitting parameters were de- 
termined by PD96 to be 



A = 0.482(1 +n/3)" 
B = 0.226(1 +n/3)" 
Q = 3.310(1 +n/3)~ 
/3 = 0.862(1 + n/3)" 
V = 11.55(1 +n/3)" 



(19) 
(20) 
(21) 
(22) 
(23) 



which are all dependent on the local slope of the linear power 
spectrum, 

dlnP , , 

-^in^- (24) 

In addition to a change of amplitude there must also be 
a change of scale. As linear perturbations become nonlinear, 
objects with density contrast 1 -f 5 have collapsed radially 
by a factor of (1 -I- (5)^' '^. This suggests that the scale /cnl is 
related to the scale /cl by the equation 



-[1- 



Xk^)] 



1/3 . 



(25) 



Taken together, Eqns. ( |17| ) and ( |25[ ) and their attending 
equations comprise the nonlinearization process, which pre- 
dicts the observed nonlinear power spectrum in redshift 
space from a given linear spectrum in real space. 

A more useful applications of these methods would be 
the inverse process - linearization. Eq. i^M can simply be 
solved for A^ (fc) and the fitting formula in Eq. (nS) can be 
inverted numerically. The change of scale expressed in Eq. 
([25I) can be stated equivalently as 



fcL = [l-f A?,L(fc^ 



-1/3 



fcNL, 



(26) 



completing the tool set for use in manipulating the power 
spectrum. 
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This method works well for scale-free power law spec- 
tra with a constant slope n of the linear P{k), but not for 
scale-dependent power spectra. The curving shape of CDM- 
type power spectra could not be reconstructed with much 
success. Jain et al. (1995) suggested that CDM-type spec- 
tra could be treated by choosing a constant 'effective' slope, 
rioff, to take the place of n in Eqns. (|l9l-C3). The effective 
slope is defined as the tangent slope of the linear P{k) at 
the nonlinear scale (where the mass fluctuation is unity). A 
linearization using ricft can only approximate the true shape 
of the linear spectrum since the method is based on deriv- 
ing the linear power from the nonlinear power on a smaller 
scale. Because the slope of the CDM-type spectrum is always 
changing, the effective index should continually change with 
scale, as was recognized by PD96. We choose n from the 
linear power spectrum on the appropriate linear scale /cl by 
numerically differentiating P{ki^), which works reasonably 
well. 



5 APPLICATION OF THE METHOD TO 
SIMULATIONS 

In order to thoroughly test this method, we applied it to 
our database of simulations or ran new simulations suit- 
able for this purpose. For each realization listed in Table H, 
we calculated the nonlinear power spectrum at the moment 
z — and then averaged them together for each box size. 
For the TACDM0.4 model, which has multiple box sizes, we 
constructed a complete nonlinear power spectrum by com- 
bining the spectra starting with the power from the largest 
box size and adding power from successively smaller boxes 
multiplied by a small factor if needed (to account for power 
from missing long waves). 

With this information we could test the procedure in ei- 
ther of two ways: linearizing the nonlinear power spectrum 
and then comparing it to the linear power spectrum of the 
model, or nonlinearizing the generated linear power spec- 
trum and then comparing it to the actual nonlinear spec- 
trum from simulations. In fact, we did both so that we could 
demonstrate the consistency of both approaches. We tested 
the linear-nonlinear corrections of §M.2 on the final nonlin- 
ear and linear power spectra of the models. The results are 
presented in Figure H for all of the models that we used 
with the exception of ACDM0.4 (since it is very similar to 
TACDM0.4). These plots incorporate the shifted mapping 
n-L(fcL) — > nL(fcL/2) from PD96. All models appear to fit 
quite well over the entire range in scales that we have been 
able to faithfully simulate. 

Figure presents the relation between A^l and Al- 
This format removes k dependence by plotting the linear- 
nonlinear mapping only. This allows a more careful exam- 
ination of the goodness of fit than the previous Figure. In 
addition to plotting the shifted linear-nonlinear mapping 
of PD96, the mapping with no shift in nL(feL) is plotted 
for comparison. We find that most models favor the shifted 
mapping except for ACDM0.4 model which clearly favored 
no shift, for which we have no explanation. 



6 APPLICATION OF THE METHOD TO 
OBSERVATIONS 

Now that a reliable linearization method has been designed 
and tested, it can be applied to real data in an attempt 
to constrain the models that best represent the Universe. 
As opposed to linearizing the observational data based on 
an arbitrary model and then comparing the linear power of 
a model in question with the linearized spectrum (PD94, 
Ma 1996), a particular model was applied to the observa- 
tional data and the linearization was performed according to 
the specific theoretical predictions of that model. The most 
correct model should be identifiable by self-consistency: an 
exceptionally good fit of the linearized observations to the 
theoretical linear spectrum. 

We adopted the same compilation of observational 
power spectra and error estimates as in PD94. For each data 
set we applied the linearization procedure for each cosmo- 
logical model to correct for redshift distortions (except in 
the case of APM data), bias, and nonlinear effects. Once 
all data sets were linearized, the collective sample of power 
spectrum data points were averaged in bins 0.1 logjQ k wide. 
The bin widths and centers were chosen to be nearly iden- 
tical to those of PD94 in their Table 1. The resulting plots 
appear in Figure |5| and are discussed below. 

Our approach in treating errors is different from PD94. 
They simply accepted errors based on counting statistics in 
bins, thus neglecting internal errors of the data sets. We at- 
tempted to keep track of individual error bars throughout 
the linearization process, assuming that errors of points in 
each data set are independent. This is not actually true, but 
it is difficult to estimate how the data points correlate. The 
error estimates were made in the following way. When per- 
forming redshift and bias corrections, the data points were 
simply scaled up or down by a factor (Eq. ([13|)), so the er- 
ror estimates were simply scaled by the same factor. For 
the case of nonlinear corrections, we performed the correc- 
tion on both A^ and the upper limit A^ + 5 and used the 
difference between these two corrected points as the error es- 
timate. Technically, this gives an error at a different k than 
the data point in question, which has the effect of reducing 
the size of the error. However, since the original errors are 
uncertain to begin with, we expect that our final estimates 
are at least a fair representation of the errors. Correlations in 
the observational data points will lead to the plotted errors 
being an underestimate of the true errors. 

While we assumed the same relative biases as PD94 
(Eq. (u)), the biases were all scaled by the same constant 
so that the linearized spectrum obtained from observational 
data could be best fit to the linear dark matter spectrum ac- 
cording to a x^ test. Each bias in Eq. (pi) was multiplied by 
this factor, reflected in the value of fei, before being applied 
to the observational data through Eq. (|l3[). Table K gives 
the best 61 for each model as well as an indicator meant to 
quantify the goodness of fit of the linearizations. The indica- 
tor (called the 'P ratio') is the average of the ratio of linear 
power to linearized power at three points near the middle of 
the sample, k =0.071, 0.090, and 0.113 h Mpc"\ chosen to 
be near the maximum curvature of the power spectra and 
have relatively small error bars. Since decreasing 61 moves 
the spectrum up and increases its curvature, this indicator 
should be a reliable estimate of how close the curvature of 
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0.01 



0.01 



0.1 1 

k (hMpc"') 



10 



0.1 1 

k (hMpc"') 



10 



Figure 3. This figure plots tlie linear power spectrum for each model (dotted fine) and the nonlinear power spectrum from simulations 
(solid line) against the linearized power spectrum (short dashed line) and the nonlinearized power spectrum (long dashed line). The 
linearized and nonlinearized lines are from the application of the method described in §W, with a modification of the local slope of the 

linear power spectrum nLC^^L) ' "lC^l/S) as used in PD96. Because the TACDM0.4 model has simulations with high resolution that 

sample smaller waves, it is plotted to higher values of k. The short horizontal lines denote boundaries between the 300, 100, and 50 
h~ Mpc samples. At high values of k resolution effects in the simulations bring the power spectrum lower than what it should be. We 
chose to stop plotting curves once A^(fc) turned over. 



the linearized spectrum is to matching that of the linear 
power spectrum for the model. In the case where the cur- 
vature is very different, however, the P ratio tends to be 
low (since the two best-fitting curves intersect at the mid- 
point) giving a false indication of a good fit. In our sample of 
models, the P ratio for the SCDM model is the only one to 



suffer from this problem because the SCDM power spectrum 
simply has the wrong shape. 

Note that the last few (i.e., highest k) binned points are 
always less reliable than the error bars indicate. They are 
either averages of only the APM and CfA results or of APM 
alone and have few points per bin. Quite often (e.g., Liddle 
et al. 1996a) the highest-fc points are ignored because the 
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Figure 4. The effects of using sliifted or unsliifted n on tiie linear-nonlinear mapping. The solid line represents the mapping in Eq. ( |l8| ) 
performed without a shift while the dotted line represents the mapping with a shift of a factor of 2, as in PD96. The filled circles show 
the linear— nonlinear relation of the real-space power spectrum from simulations of the models. 



earlier linearization algorithms were not accurate for these 
points. However, we think that these points should be in- 
cluded in the analysis as long as there are enough data points 
in the bins. They have the smallest statistical error bars and 
now cannot be excluded solely on the basis of failure of the 
linearization procedure, since our method demonstratably 
works up to this range in k. 

We find that the ACDM family of models have the 
wrong shape when the best fitting bias bi is determined. 
All of these models have too much curvature and fit poorly 
at the high-fc tail of points, even though most of the points 



are within the nominal error bars. As bi is lowered below 
unity, the curves of the reconstructed linear power spectra 
rise in amplitude and increase in curvature, which is what 
is responsible for the relatively poor fits. The reconstructed 
CHDM models both seem to do an exceptionally good job of 
fitting their respective power spectra. The excellent match 
of the curvature and the small adjustment to the bias make 
CHDM the most successful class of candidate models in our 
sample. 

Because the ACDM power spectra did not fit well, we 
tried examining the effects of another possible modification 
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Figure 5. The application of the linearization technique to the observational data set. The linear power spectrum of the imposed model 
(dotted line) and two different cases of bias parameters are plotted (dots with error bars). The open dots correspond to the bias level of 
PD94 (Eq. (h) in this paper) and the filled dots are for the bias level adjusted by the factor given in Table bl as determined by minimizing 
X^ ■ Error bars represent internal errors of each point in the bin added in quadrature. See the text for a more detailed explanation. 



process. The formalism of Ml assumes that the bias does 
not change with scale. This is a useful first approximation, 
but may prove to be an oversimplification. It is quite possi- 
ble that the bias depends on scale. Preliminary results from 
sophisticated A''-body plus hydro simulations of galaxy for- 
mation done by Yepes et al. (1996) (which include multi- 
phase treatment of gas and supernovae feedback) indicate 
that bias is slightly rising with k for both CDM and ACDM 
models. The semi-analytic merging hierarchy inclusion of 
hydrodynamic effect in simulations by Kauffmann, Nusser, 



& Steinmetz (1995) gave similar results regarding scale de- 
pendence of bias. In order to obtain a rough estimate of 
the possible effect we introduce a very simple model for bias 
which modifies the redshift-space correction of Eq. (113) . Fig- 
ure H shows the form of a bias which increases with scale in 
(a) and which decreases with scale in (b). It is a linear func- 
tion of fc: bi = 1 -I- fcNL/(l h Mpc~^) which is a factor of two 
at fcNL = 1 h Mpc~^. Panel (c) shows that an extra positive 
bias makes the best fit even worse while (d) shows that anti- 
bias improves the best fit by raising the high-fc tail. Because 
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Table 3. Summary of best-fit biases 
Model 



b^ 



l.SfeJ 



P ratio'* 



SCDM 


1.42 


1.85 


1.35 ±0.24 


0.98 


ACDMo.3 


0.50 


0.64 


0.69 ±0.12 


1.41 


ACDMo.4 


0.58 


0.75 


0.81 ±0.14 


1.34 


TACDMo.4 


0.69 


0.90 


0.81 ±0.14 


1.33 


ACDMo.5 


0.85 


1.11 


0.92 ±0.16 


1.09 


CHDMo.7 


0.87 


1.13 


_ 


0.99 


CHDMo.8 


0.87 


1.13 


- 


1.04 



"• Factor by which the PD94 bias must be multiplyed to achieve 
best fit. 

Bias of blue galaxies as indicated by the best fit fej. 
'^ Bias of blue galaxies from POTENT analysis Hudson et al. (1995). 
'^ Average of PiinearZ-Plincarizod at k =0.071, 0.090, and 0.113 
h Mpc-i. 




0.1 
k (ftMpc"') 



Figure 6. Test forms of scale-dependent bias and their effects on 
the recovered linear power spectrum for TACDM0.4. The solid 
line is the redshift correction of Eq. (t3), the dotted line is the 
scale-dependent bias function and the dashed line is the applica- 
tion of the function to the redshift correction, (a) Bias increasing 
linearly with scale to a factor of 2 at fc = 1. (b) Bias decreasing 
linearly with scale to a factor of 2 at fe = 1. (c) The positive bias 
increases the curvature requiring more antibias. (d) The antibias 
decreases the curvature requiring less antibias. 



such a large anti-bias is difficult to motivate physically, we 
do not regard it as an acceptable solution to this problem 
(cf. KPH96). 

At this point we can demonstrate that the choice of 
ID velocity dispersion has an impact on the form of a re- 
constructed linear power spectrum, as was ari 



in 5 4.1 



Figure ^ shows that the velocity dispersion separates the dif- 
ferent cases at the high-fc tail for the unbiased TACDM0.4 
model. Three realistic values of a were chosen, 300 km s~^ 
as in PD94, 430 km s~^ as was used for all ACDM models, 
and 560 km s~^ which is representative of CHDM models. 
Since for k > 0.2 h Mpc~^ the curves are separated by an 



TACDM„, 
b = 1.0 

1 = 300 km/s 

cr=430 km/s 

CT = 560 km/s 




k (iiMpc"') 

Figure 7. The effect of different choices of ID velocity dispersion 
on the reconstructed observational TACDM0.4 power spectrum. 
Error bars have been included for the a = 430 km s~^ case only, 
but are representative of the other cases. Only the high-fc points 
are shown since the differences are minimal at low k (see Figure ih. 
The last two points are degenerate for all values of a since they 
are constructed from APM data only, which are not affected by 
the redshift correction. 



amount greater than the error bars, the redshift correction 
can have a noticeable effect on the amplitude of the recon- 
structed power spectrum. 

An interesting thing that we noticed was that the recon- 
structed observational power spectra of our ACDM models 
were virtually identical. The TACDM0.4 model was the only 
exception (because of the tilt), but the difference is small. 
Figure n, shows the average of the ACDM0.3, ACDM0.4 and 
ACDM0.5 models along with the TACDM0.4 model for 61 
values 1.0, 0.8 and 0.6. Using this to analyze the fit of 
the general ACDM spectrum, the averaged reconstructed 
ACDM spectrum for the same three bi values overlaid by the 
linear power spectra of all of the ACDM models are shown 
in Figure H. Again, this figure demonstrates that no level of 
constant 61 can provide an adequate reconstruction of the 
linear power spectrum for the ACDM models: the shape of 
the spectrum is not optimal. 

A similar finding is that there is a universal linear power 
spectrum reconstructed from the observational data. Re- 
construction for all models that we examined resulted in 
practically the same spectrum in the range 0.01 ^ k ^ 0.1 
h Mpc^^. Figure hd shows that all models share an almost 
identical shape and amplitude in this region. This means 
that any future progress in distinguishing between these 
models will need to be emphasize nonlinear scales k <; 0.2 
h Mpc~^. The data shown in Figure hd have been tabulated 
in Tables |-|[ 

They contain the linearized observational power spectra 
for the ACDM, CHDM, and SCDM models, respectively. 
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Table 4. ACDM family of recovered linear power spectra A^(fc) : 



,P{k) 





fei 


= 0.6 


h 


= 0.8 


bi 


= 1.0 


fc,(fcMpc-i) 


ACDM 


TACDM 


ACDM 


TACDM 


ACDM 


TACDM 


0.0142 


0.0028 


0.0028 


0.0016 


0.0016 


0.0010 


0.0010 


0.0179 


0.0033 


0.0033 


0.0020 


0.0020 


0.0013 


0.0013 


0.0225 


0.0082 


0.0082 


0.0050 


0.0049 


0.0033 


0.0033 


0.0284 


0.0234 


0.0234 


0.0145 


0.0145 


0.0098 


0.0098 


0.0357 


0.0641 


0.0638 


0.0400 


0.0400 


0.0274 


0.0273 


0.0450 


0.0935 


0.0932 


0.0583 


0.0582 


0.0397 


0.0397 


0.0566 


0.156 


0.155 


0.0970 


0.0969 


0.0597 


0.0597 


0.0713 


0.275 


0.269 


0.179 


0.178 


0.121 


0.121 


0.0897 


0.461 


0.456 


0.270 


0.273 


0.182 


0.182 


0.113 


0.639 


0.624 


0.387 


0.387 


0.248 


0.249 


0.142 


0.942 


0.921 


0.621 


0.611 


0.406 


0.403 


0.179 


1.17 


1.14 


0.857 


0.845 


0.593 


0.615 


0.225 


1.72 


1.63 


1.15 


1.16 


0.821 


0.792 


0.284 


1.53 


1.33 


1.39 


1.30 


1.17 


1.16 


0.357 


1.80 


1.62 


1.37 


1.22 


1.36 


1.25 


0.450 






1.67 


1.52 


1.26 


1.17 


0.566 










1.49 


1.36 



Table 5. CHDM recovered linear power spectra A'^(fc) : 



,P{k) 





h = 


:0.8 


bi = 


: 1.0 


k (h Mpc-1) 


CHDMo.7 


CHDMo.8 


CHDMo.7 


CHDMo.8 


0.0142 


0.0016 


0.0016 


0.0010 


0.0010 


0.0179 


0.0018 


0.0018 


0.0012 


0.0012 


0.0225 


0.0044 


0.0044 


0.0030 


0.0030 


0.0284 


0.0122 


0.0122 


0.0085 


0.0085 


0.0357 


0.0336 


0.0336 


0.0237 


0.0237 


0.0450 


0.0500 


0.0499 


0.0350 


0.0350 


0.0566 


0.0863 


0.0861 


0.0530 


0.0529 


0.0713 


0.159 


0.158 


0.111 


0.111 


0.0897 


0.236 


0.236 


0.170 


0.169 


0.113 


0.341 


0.350 


0.235 


0.233 


0.142 


0.537 


0.615 


0.368 


0.388 


0.179 


0.647 


0.813 


0.540 


0.619 


0.225 


0.768 


1.11 


0.666 


0.854 


0.284 


0.797 


1.15 


0.768 


1.11 


0.357 


0.762 


1.06 


0.690 


0.868 


0.450 


0.830 


1.26 


0.750 


1.03 


0.566 






0.795 


1.16 



7 CONCLUSIONS 

We have reconstructed the linear power spectra from ob- 
servations for several current cosmological models based on 
rigorously tested linearization methods. We have found that: 

(i) The linearized power spectra of all ACDM and 
CHDM models are systematically below their respective lin- 
ear power spectra using the bias ratio in Eq. (pi) normalized 
so that the bias of IRAS galaxies bi = 1.0. It was neces- 
sary to adjust the normalization of the bias ratio so that 
fei < 1.0 to obtain the best fit. While this implies that op- 
tically selected galaxies must be antibiased (6o < 1) to be 
consistent with models, this may be acceptable since a re- 
cent POTENT analysis (Hudson et al. 1995) allows a slight 
antibias for some ACDM models. 

(ii) The linearization of the SCDM model cannot be 
made to fit its linear power spectrum under any circum- 
stances and fits very poorly at best. This simply reconfirms 
the incompatibility of SCDM with observations, as has long 
been recognized. 



(iii) The ACDM models all appear moderately success- 
ful, as indicated by our 'P ratio' indicator in Table H, but 
share difficulties in matching the shape of their linear spec- 
tra. This is a consequence of our finding that the linearized 
spectra of all ACDM models share a common shape for a 
given bias. The only exception is the TACDM0.4 model, 
though it deviates from the others by at most 10 per cent at 
fc ~ 0.5 ft Mpc"^ This is to be expected since TACDM0.4 
has more curvature by design. The best fitting ACDM model 
is ACDM0.5 as determined from Table H. 

(iv) The CHDM models seem to do an excellent job 
of recovering their linear power spectra from observations. 
They require relatively small 61 to fit well and have the low- 
est P ratios. Ma (1996) found that a slight tilt of n = 0.9- 
0.95 is needed to bring h = 0.5 CHDM models into agree- 
ment with the PD94 reconstructed spectrum. However, if 
one reconstructs the linear spectrum self-consistently, using 
the nonlinear spectrum for the CHDM models considered 
here, there is good agreement without any tilt. 

(v) There are physical motivations to include a scale- 
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Figure 8. Observational linear power spectrum reconstructed 
under an assumption that the Universe is of a ACDM type. Be- 
cause the reconstructed spectra were very similar for ACDM0.3, 
ACDM0.4, and ACDM0.5 models, we present averaged spectra for 
assumed level of bias. The spectrum reconstructed assuming the 
TACDM0.4 models differ slightly from the rest and are drawn as 
the dashed lines. These data points are given in Table U. 



Table 6. SCDM 



recovered linear power spectra A^(A:) 



k {h Mpc-i) fei 



0.6 



0.8 



1.0 



0.0142 


0.0028 


0.0016 


0.0010 


0.0179 


0.0028 


0.0018 


0.0012 


0.0225 


0.0071 


0.0044 


0.0030 


0.0284 


0.0189 


0.0122 


0.0085 


0.0357 


0.0513 


0.0335 


0.0236 


0.0450 


0.0765 


0.0498 


0.0349 


0.0566 


0.135 


0.0861 


0.0527 
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Figure 9. Comparison of averaged ACDM0.3, ACDM0.4 and 
ACDM0.5 simulations with linearized observational power spec- 
tra appropriate to ACDM cosmology for different biases. Dotted 
lines depict the linear power spectrum of ACDM0.3, short dashed 
lines depict ACDM0.4, long dashed lines depict TACDM0.4, and 
dot-dashed lines depict ACDM0.5. The biases in the panels are 
(a) bi = 1.0, (b) bi = 0.8, and (c) 61 = 0.6. 



dependent bias in the treatment of observational power spec- 
tra. Such a bias would increase with higher k and would 
make the reconstructions at fe ^ 0.2 h Mpc~^ worse than 
they already are. A bias that decreased with scale could sig- 
nificantly improve the agreement between the assumed and 
reconstructed power spectrum for each model, especially the 
ACDM models, but this appears to be an entirely ad hoc as- 
sumption. 

(vi) There is a unique model-independent linear power 
spectrum that can be recovered from observations in the 



region 0.01 Ss k ^ 0.1 h Mpc ^, at least for the models 
considered here. 

We have demonstrated that the linearization method is 
accurate to within 20 per cent up to linear scales of at least 
k = 1 h Mpc~^ and possibly higher. Because the most favor- 
able models are nearly indentical up to the highest reliable 
bin, a power spectrum with small errors that extends from 
large nonlinear scales to fc ~ 10 /i Mpc~^ or more would be 
very helpful in constraining the shape of theoretical linear 
power spectra. 
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power spectra with h\ = 1.0 in each case. 
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APPENDIX A: FITTING FUNCTIONS FOR 
LINEAR POWER SPECTRA 

Al Calculations 

The linear calculations were done using the techniques of 
Holtzman (1989). Perturbations in radiation, baryons, cold 
dark matter, and massive and massless neutrinos were calu- 
clated from a scale factor of 1 x 10~^° to the present (where 
the scale factor is a = 1/(1 + z)). For models without mas- 
sive neutrinos, three massless neutrino species were used; for 
models with massive neutrinos, the mass was either placed 
in a single neutrino species or two equal-mass species, with 
the other species remaining massless. 

For the radiation, the angular dependence of the per- 
turbations was handled by expanding the perturbation in 
Legendre polynomials, &^ = ^ Ai(fc, a)P;(cos0). The cal- 
culation started with a small number of orders, but addi- 
tional orders were added as they became needed to keep the 
final P(fc) results accurate to 1 per cent. For large scales 
(fc < 0.7 h Mpc~^), the radiation perturbations were inte- 
grated all the way to the present. For smaller scales, they 
were integrated through recombination, and then the ana- 
lytic solution presented by Bond and Efstathiou (1987) was 
used to calculate the perturbations at the present time; no 
approximations were made for any scale for the radiation 
perturbations. The recombination history of Peebles (1967) 



was used for hydrogen, but a helium mass fraction of 0.25 
was included. 

For the massive neutrinos, the perturbations are a func- 
tion of both angle and neutrino momentum. We calculated 
the evolution for 15 separate momenta up to the time when 
the neutrinos became nonrelativistic; this was done using the 
integro-differential equations presented by Bond and Szalay 
(1983). This technique allows only the 0th and 1st orders of 
the neutrino angular distribution to be computed without 
following higher orders; these are the only orders than enter 
as the graviational source terms. For small scales, even this 
computation can get very expensive, so the massive neu- 
trino perturbations were manually damped when they be- 
came gravitationally negligible (which occurs because the 
perturbations are destroyed by free-streaming). When the 
neutrinos became nonrelativistic, the full angular distribu- 
tion for each momenta was calculated, and the subsequent 
evolution was computed using the full set of coupled dif- 
ferential equations exactly analagous to those used for the 
radiation perturbation. 



A2 Fitting functions 

For convenience, we provide fits for the power spectrum of 
the hot and cold components as a function of wavenumber, 
fc. For all of the calculations presented in this paper, we used 
the computed spectra directly, rather than the fits, to obtain 
results. The fits are good to better than 5 per cent in total 
power, in the worst case. Errors of 1-2 per cent are more 
typical for k — 0.1-30 h Mpc~^. For fitting functions, we 
used 

Pc+b(fc) 



(1 + aifci/2 + a2k + agP/^ 4- a^k'^f 
for the cold and baryonic matter, and 

p .^^ ^ Pc+b(fc)exp(-cfc^/^) 

(1 -I- difcV2 + dafc + dsfc^/^ 4^ d4fc2)^ 

for the massive neutrinos. The total power is given by 



(Al) 



(A2) 



P(fc) = 



Sib + flc 

fio 



"0 



"Kik) 



(A3) 



The parameters for all the models considered in this 
paper are presented in Table [H for the cold component and 
in Table p for the hot component; the coefficients ass ume 



the use of units of h Mpc for k in equations Al 



and 



A2 



The normalization A is set to reproduce the cts used in this 
paper. 
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